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Abstract 

Recent studies exploring the nnolecular genetic basis for nnigratory variation in 
aninnals have identified polynnorphisnns in two genes (CLOCK and ADCYAP1) 
that are linked to circadian rhythnns and correlate with nnigratory propensity and 
phenology annong individuals and populations. Results fronn these initial studies 
are nnixed, however, and additional data are needed to assess the generality 
and diversity of the nnolecular nnechanisnns that regulate the biology of 
nnigration. We sequenced CLOCK and ADCYAP1 in 15 populations across the 
two species of the avian genus Junco, a North Annerican lineage in which 
nnultiple recently diverged subspecies and populations range fronn sedentary to 
long-distance nnigrants. We found no consistent associations between allele 
length and nnigratory status across the genus for either CLOCK or ADCYAP1 . 
However, within two subspecies groups, populations that nnigrate longer 
distances have longer CLOCK alleles on average. Additionally, there was a 
positive relationship between ADCYAP1 allele length and nnigratory 
restlessness (zugunruhe) annong individuals within one of two captive 
populations studied — a result sinnilar to those reported previously within captive 
blackcaps (Sylvia atricapilla). We conclude that, while both ADCYAP1 and 
CLOCK nnay correlate with nnigratory propensity within or annong certain 
populations or species, previously identified relationships between nnigratory 
behavior and sequence variants cannot be easily generalized across taxa. 
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Introduction 

Every year billions of birds make round-trip flights from their 
breeding grounds to suitable winter climes and back again ^ This 
seasonal migration requires immense coordination of different sys- 
tems including fattening, locomotion, orientation, and activity-. 
Several aspects of migratory behavior including onset, duration, 
intensity, and orientation are known to be heritable and to respond 
quickly to artificial and natural selection - indicating a degree of 
genetic control^"^. However, other studies indicate that both learn- 
ing and developmental environments can also shape migratory 
phenotypes-'^'*^. 

The propensity to migrate appears extremely labile on an evolu- 
tionary scale, with shifts between migratory and sedentary status 
independently arising, often repeatedly, in multiple lineages'""^-. In 
recent decades, many populations and species are reducing the dis- 
tance or frequency of migrations, and others are ceasing to migrate 
altogether, presumably in response to changing environments 
world-wide^' ^^"^^ Understanding the molecular genetic mecha- 
nisms that regulate migration biology is imperative both for 1) 
understanding the evolutionary processes of behavioral adaptation 
and diversification, and 2) conserving and managing the phenom- 
enon of animal migrations in the face of ongoing environmen- 
tal change. However, relatively few studies have addressed the 
molecular genetic mechanisms of migration^^. Our goal in this 
study was to build upon previous candidate gene studies to deter- 
mine whether or not similar genetic mechanisms explain variation 
in migratory propensity at the level of individuals, populations, and 
species across taxa. 

Previous studies have identified a number of candidate genes that 
may be involved in regulating migration, with particular focus on 
circadian rhythm genes. For example, the gene CLOCK controls 
several aspects of the circadian rhythm in mice^^, and in the blue 
tit (Cy artistes caeruleus) it varies across latitudinal clines^^ and 
with timing of breeding In salmon (Oncorhynchus spp.), CLOCK 
microsatellite repeat length also varies geographically and predicts 
migratory timing-^^"^^. Monarch butterflies (Danaus plexippus) use 
CLOCK as part of the circadian pathway that guides their annual 
migration-^ suggesting that CLOCK may play a role in the timing 
of migration across vertebrate and invertebrate lineages. However, 
blackcaps (Sylvia atricapilla^^) bluethroats {Luscinia svecica^^), 
and swallows of the genus Tachycineta-^ do not show such correla- 
tions between this CLOCK polymorphism and migratory behavior 
or geographic location, despite possessing genetic variation at the 
locus. Furthermore, barn swallows (Hirundo rustica) show substan- 
tial variation in migratory phenotype, but very little variation at this 
CLOCK polymorphism-^'. 

In another gene, adenylate cyclase-activating polypeptide 1 
{ADCYAPl), length of a microsatellite repeat predicts migratory pro- 
pensity of both populations and individuals in blackcaps-^. ADCYAPl 
is expressed throughout the brain and body of vertebrates-^ and 
encodes pituitary adenylate cyclase-activating polypeptide (PACAP) 
(reviewed in-^). PACAP is involved in many diverse behavioral and 
physiological phenotypes, including the circadian system (reviewed 
by Vaudry et al}^). Specifically, in chickens {Callus gallus), PACAP 
directly activates CLOCK and other circadian genes in the pineal 
gland-'^. Further, PACAP stimulates release of melatonin^", and is a 



key component in entraining the circadian system to the light-dark 
cycle^^ suggesting that it may be involved in circannual rhythms that 
rely on sensing day length, such as migration. Additionally, PACAP 
is involved in a large number of physiological and behavioral effects, 
from feeding behavior to breathing^^ many of which may be related 
to physiological changes involved in migration-^. 

The body of work to date thus indicates that certain genes are 
associated with migratory phenotype in some way across phyla and 
within classes. However, these findings also warrant further study 
of the genetics of this complex phenotype at multiple levels in inde- 
pendent systems to determine whether selection is using conserved 
genetic machinery to arrive at similar phenotypic outcomes (e.g., 
Korsten et aV^-). 

Here, we examined allelic variation in both CLOCK and ADCYAPl 
in the genus Junco at the levels of species, populations, and individu- 
als. The junco species group is currently comprised of the yellow- 
eyed junco (/. phaeonotus) and the dark-eyed junco (/. hyemalis), 
each of which contains multiple subspecies and populations showing 
a wide range of migratory behaviors^^"^^ Several of these subspecies 
have been identified as distinct species in the past (e.g. Miller^^) and 
investigations into the genetic structure and systematics of the group 
are ongoing (e.g. Mila et aV^\ McCormack et al?\ Rasner et aV^ 
and Whittaker et aV"^). 

The various dark-eyed junco subspecies offer an especially exciting 
opportunity to investigate migratory differences between closely re- 
lated populations. There are 15 subspecies comprising the species^"^, 
and genetic similarity indicates that these closely related forms may 
have radiated just within the last 10,000 years-spreading across 
North America following the most recent glacial maximum^^\ This 
rapid diversification has limited the ability to identify clear phyloge- 
netic patterns within the species^^'. These subspecies and populations 
differ markedly in migratory phenotype, ranging from completely 
sedentary to those that migrate thousands of kilometers (Table 1), 
with several populations also exhibiting individual variation in 
migratory behavior in which some individuals migrate and others 
do not^-'^^'^^ '^^^ Further, a coastal population of the Oregon junco 
(/. h. thurberi) has diverged from a montane altitudinal migrant popu- 
lation in the last 30 years and become sedentary in a milder environ- 
ment on the campus of the University of California in San Diego^^'^^. 

The diversity of migratory behavior along with the close genetic 
relationships among divergent subspecies and populations make the 
junco species an excellent system in which to further examine the 
role of candidate genes linked to migration. Based on earlier studies 
of CLOCK and ADCYAPl variation^ we predicted that junco 
populations that migrate longer distances would possess longer 
microsatellite repeat-length alleles of both CLOCK and ADCYAPl 
when compared with those that migrate shorter distances including 
altitudinal migrants and sedentary populations. We also predicted 
that allele length would positively covary with individual variation 
in migratory restlessness within populations. 

Methods 

Classifying population migratory behavior 

Populations of juncos were classified according to a scale of migra- 
tory distance based on published reports in the literature^^"^^''^°"^^''^'^'^^ 
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Table 1. Sampled populations and assigned migratory scores. Group categories follow generally from distinct groups of subspecies 
or unique forms as discussed by Mila et al. (2007f^ or Nolan et al. (2002)^^. Subspecies designations follow Miller (1941)^^ as summarized 
by Nolan et al. (2002)^^ and Sullivan etal. (1999)^^; dark-eyed subspecies are noted as J. h. spp. and yellow-eyed subspecies as J. p. spp. 
Migratory Behavior codes indicate the range of migratory behaviors inferred from breeding and wintering distributions, as follows: LD2 
= long distance II, at minimum 1600 km up to 5600 km, depending on migratory connectivity; LD1 = long distance I, likely 400-700 km, 
possibly up to 5000 km, depending on migratory connectivity; R = regional, typically greater than 200 km; A = altitudinal, typically less 
than 200 km; P = apparent partial migration documented; F = apparent facultative migration documented; S = sedentary (see text for 
references). Migratory score ordinally ranks population migratory propensity, similar to previous published accounts^^. 



Group 


Subspecies 


Sampling site ^^^^^^^ 


Range 


Abbreviation 


Migratory 
behaviors 


Migratory ^ 
score 


Slate-colored 


J. h. 


hyemalis 


Mississippi, USA 


Wintering 


SCJU-MS 


LD2 


6 


Slate-colored 


J. h. 


hyemalis 


Indiana, USA 


Wintering 


SCJU-IN 


LD1 


5 


Slate-colored 


J. h. 


hyemalis 


Michigan, USA 


Wintering 


SCJU-MI 


Rto LD 


4 


Carolina 


J. h. 


carolinensis 


Virginia, USA 


Breeding 


CRJU 


A, P, F 


2 


White-winged 


J. h. 


aikeni 


South Dakota, USA 


Breeding 


WWJU 


Rto LD1; A, P, F 


3 


Oregon 


J. h. 


oreganus 


British Columbia, Canada 


Breeding 


ORJU-BC 


Rto LD1; P, F 


3 


Oregon 


J. h. 


thurberi 


Mt. Laguna, California, USA 


Breeding 


ORJU-LM 


A, P, F 


2 


Oregon 


J. h. 


thurberi 


San Diego, California, USA 


Year-round 


ORJU-SD 


S 


1 


Gray-headed 


J. h. 


caniceps 


Utah, USA 


Year-round 


GHJU 


R 


3.5 


Pink-sided 


J. h. 


mearnsi 


Wyoming, USA 


Breeding 


PSJU 


R to LD1 


4.5 


Guadalupe 


J. h. 


insularis 


Guadalupe Island, Mexico 


Year-round 


GUJU 


S 




Yellow-eyed 


J. p. 


phaeonotus 


Durango, Mexico 


Year-round 


YEJU-DO 


S 




Yellow-eyed 


J. p. 


phaeonotus 


Mexico City, Mexico 


Year-round 


YEJU-DF 


S 




Guatemala 


J. p. 


alticola 


Huehuetenango, Guatemala 


Year-round 


GTJU 


S 




Baja 


J. p. 


bairdi 


Baja California Sur, Mexico 


Year-round 


BAJU 


S 





personal observation. Each population was ranked on a scale from 1 
(sedentary) to 6 (long-distance migrant) based on the average dis- 
tance migrated by individuals in that population (Table 1). Distance 
migrated also varies within populations. For example, in the Eastem 
United States, females^^' and adults^^ generally migrate farther 
from their breeding ranges than males and yearlings, respectively. 
Some junco populations also exhibit partial migration, in which 
a subset of individuals remains on the breeding grounds while 
others depart^^'^^ The classification system employed here broadly 
categorizes the general migratory phenotype of each population 
and was organized to resemble previous reports for other species 
(e.g. blackcaps^^). 

Assessing the migratory behavior of individuals 

Individuals from two populations of /. h. thurberi residing in southem 
California were captured as recently independent offspring and held 
in captivity under identical conditions in a 'common garden'. Indi- 
viduals were held individually (0.61 m x 0.61 m x 0.61 m cages) 
in a climate controlled indoor aviary with ad libitum access to food 
and water. Individuals were neither visually nor acoustically isolated, 
and light regimes were changed bi-weekly to match photoperiod at 
their home latitude. All methods were reviewed and approved by 
the lACUC at Indiana University (protocol #09-037). To reduce the 
probability of capturing siblings and pseudoreplication, individuals 
were captured from multiple locations, and one individual from each 
identified sibship (one in each population) was randomly included in 
all analyses, as in previous studies^^. Additional studies of these cap- 
tive common garden populations are reported elsewhere (^^; Atwell 
et al. in review). 



Briefly, spring nocturnal restlessness behavior was scored from 
3 March to 1 August 2010 by recording intervals of night-time 
perch-hopping activity of birds housed indoors in individual cages 
equipped with microswitched perches. This method is considered 
effective to capture 95% of migratory behavior observed in infrared 
video and ultrasound methods"^^, and even non-traditional songbird 
migrations tend to occur at night^^^'^^ Juncos typically migrate at 
night^^, and this sampling period overlapped with the typical spring 
migration as well as breeding seasons of dark-eyed juncos^^'^^ As 
was found in prior studies of migratory restlessness in the junco 
(e.g., Ketterson and Nolan^^) and other species, migratory rest- 
lessness extended into the typical breeding period, perhaps due to 
caged birds' inability to perform reproductive activities^ \ Nocturnal 
hops were recorded starting 30 minutes after dusk until 30 minutes 
before dawn under a daylight regime that simulated the photoperiod 
at their native latitude. GraphPad Prism 3.0 (GraphPad, USA) was 
used to quantify the seasonal intensity of individual spring migra- 
tory restlessness by calculating the area under each fitted seasonal 
profile curve, which provided a total migratory propensity score for 
each bird (n = 36^^; Atwell et al. in review). 

Sampling 

Analysis of population differences in genotype was conducted 
using DNA collected between 1996-2010 from adult juncos that 
were captured from breeding or wintering populations as a part 
of other studies. Less than 100 |iL of whole blood was collected 
via capilary tube from a small puncture in the alar vein for future 
extraction. Both males and females were included, and in all 
populations in which sexes were known, the sexes did not differ 
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in allele length for either locus (t-test, all p > 0.1). Birds sam- 
pled in their winter range were captured from December through 
February, after j uncos are thought to have arrived at their winter- 
ing grounds and before they have departed in the spring^"^'^"^. All 
individuals from breeding locations were captured as adults dur- 
ing the breeding season of that population. No known siblings 
were included in these analyses. 

Genotyping 

We extracted DNA from whole blood from 15 populations spanning the 
range of migratory phenotypes in the genus Junco (Table 1). We devel- 
oped primers surrounding a microsatellite repeat of interest based 
on previously pubUshed reports for CLOCK'^^' and ADCYAPP^^' 
in birds, and we modified the primers to more closely match the 
junco sequence based on the junco transcriptome (^^; See Table 2 
for primer details). The primers were used to amplify the respective 
polymorphisms in a single multiplex reaction using an amplifica- 
tion kit and following the manufacturer's directions (Qiagen Inc, 
Valencia, California, USA). An initial 15' heat activation cycle 
(95°C) was followed by 35 cycles of 94°C denaturation for 30", 
60°C annealing for 90", and 72°C extension for 60", and a final 
elongation of 10^ PGR products were then diluted 1:200 in ddH^O 
with LIZladder (Applied Biosystems, Carlsbad, California, USA), 
denatured at 95°C for 5^ and placed on ice. Samples were analyzed 
on an ABI 3730 using Peak Scanner vl.O (Applied Biosystems, 
Carlsbad, California, USA). 

Data analysis 

In order to assess the relationship between allele length and 
migratory behavior among populations for both loci, we performed 



Spearman's correlations between population mean allele length 
and population migratory status. Qualitatively similar results are 
obtained with major allele count methods that have been applied in 
previous studies ^^'^^ and when only dark-eyed subspecies were con- 
sidered (data not shown). Within the sub-species groups for which 
we sampled multiple populations (Oregon juncos and slate-colored 
juncos, see Table 1), we additionally performed a nested ANOVA 
on allele lengths with population as a grouping factor within each 
subspecies to determine whether mean allele length differed signifi- 
cantly between populations or these subspecies. 

To relate individual variation in allele length of ADCYAPl and 
CLOCK to migratory restlessness, we fit linear models with migra- 
tory restlessness as the dependent variable and all combinations of 
population, mean allele length of the individual, sex, and each pair- 
wise interaction as predictors. The resulting models were compared 
against each other using Akaike's Information Criterion. When the 
best fit model yielded a significant interaction term, we followed 
with a separate linear model for each populations^'. 

All statistical analyses were performed in R version 2.15.2^^. 
Hardy-Weinberg equilibrium and heterozygosity were calculated 
with the package genetics^l All reported p-values are two-tailed, and 
hence conservative with respect to a priori directional hypotheses. 

Results 

Among population comparisons for ADCYAP1 

We identified 1 6 alleles at the ADCYAPl locus ranging in length from 
154 to 181 base pairs (bp), and all populations contained between 
five and ten alleles (Table 3). Allele lengths of 161 or 163 bp were 



Table 2. Primers used in this study. Both primers are based on Steinmeyer et al^^, and CLOCK was 
modified to matcli junco sequences^^. 



Gene name 


Forward primer 5' 3' VHHI 


W§ Reverse primer 5' 3' SHUI^^I 


CLOCK 


TTTTCTCAAGGTCAGCAACTTGT 


CTGTAGGAACTGCTGGGGKTGCTG 



ADCYAP1 GATGTGAGTAACCAGCCACT ATAACACAGGAGCGGTGA 



Table 3. Distribution ADCYAP1 allele length within each population. Table lists the number of individuals from each population 
possessing each of the identified ADCYAP1 alleles. Number of individuals genotyped is indicated as "n" for each population. 
Het = calculated heterozygosity for the population. Other abbreviations follow from Table 1 . 



Population 


154 


155 


156 


157 


158 


159 


160 


161 


162 


163 


164 


165 


167 


169 


171 


181 


n 


Het 1 


SCJU-IN 


1 


3 


0 


3 


0 


9 


0 


11 


0 


8 


0 


8 


1 


1 


1 


0 


19 


0.809 


SCJU-MI 


0 


1 


0 


1 


0 


11 


1 


11 


2 


17 


0 


11 


2 


1 


0 


0 


23 


0.852 


SCJU-MS 


0 


0 


0 


4 


0 


6 


0 


8 


0 


13 


0 


4 


2 


1 


0 


0 


29 


0.816 


CRJU 


0 


0 


0 


2 


0 


9 


0 


19 


0 


12 


0 


1 


1 


0 


0 


0 


22 


0.710 


WWJU 


0 


0 


1 


5 


0 


1 


0 


35 


0 


7 


0 


5 


1 


1 


0 


0 


28 


0.587 


ORJU-BC 


0 


1 


0 


1 


0 


3 


0 


12 


0 


8 


0 


2 


3 


0 


0 


0 


15 


0.768 


ORJU-LM 


0 


0 


0 


1 


0 


11 


0 


18 


0 


16 


5 


14 


1 


0 


0 


0 


32 


0.802 


ORJU-SD 


0 


0 


0 


10 


0 


8 


0 


26 


0 


15 


0 


3 


12 


0 


0 


0 


36 


0.791 


GHJU 


0 


0 


0 


2 


0 


5 


1 


17 


0 


17 


0 


8 


4 


0 


0 


0 


27 


0.778 


PSJU 


0 


1 


0 


4 


0 


2 


0 


16 


0 


14 


0 


7 


2 


2 


0 


0 


24 


0.786 


GUJU 


0 


0 


0 


0 


0 


6 


0 


11 


0 


6 


0 


8 


4 


0 


0 


1 


18 


0.811 


YEJU-DO 


0 


3 


0 


2 


0 


4 


0 


10 


0 


1 


0 


0 


6 


2 


0 


0 


14 


0.812 


YEJU-DF 


0 


0 


0 


2 


1 


18 


1 


17 


0 


13 


0 


4 


6 


0 


0 


0 


31 


0.794 


GTJU 


0 


0 


0 


4 


0 


8 


0 


13 


0 


0 


0 


7 


2 


0 


0 


0 


17 


0.761 


BAJU 


0 


0 


0 


0 


0 


2 


0 


3 


0 


21 


0 


16 


3 


1 


0 


2 


24 


0.700 
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most common in all populations. Nearly all of the identified alleles 
differed by multiples of two bases from these common alleles. The 
exceptions were represented in only 12 individuals, including five 
individuals from the Laguna Mountain population that shared an 
allele 164 bp long. All populations, both separately and combined, 
were in Hardy-Weinberg equilibrium (all p > 0.73). 

There was no correlation between the population mean allele 
length of ADCYAPl and population migratory status (Spearman's 



correlation, rho = 0.022; p = 0.938; Figure 1). Within the sub-species 
groups sampled at multiple locations the length of the ADCYAPl 
allele did not differ between subspecies, or within subspecies (nest- 
ed ANOVA, subspecies difference: - 0-233, p = 0.629; popula- 
tion within subspecies: ^^296 ~ 0-852, p = 0.493; Figure 2). 

Among population comparisons for CLOCK 

We identified eight alleles at the CLOCK locus that varied in length 
from 267 to 285 bp, and all but one population (/. h. insularis 



Q_ 

O 
Q 
< 



TO 
3 
Q. 
O 
Q_ 



BAJU 
O 



GUJU 
O 



ORJU-LM 
O 



GHJU 
O 



SCJU-MI 



YEJU-DO 
O 

YEJU-DF 
O 
GTJU 



CRJU 
O 



O 

ORJU-BC 
WWJU 
O 



SCJU-MS 
O 



SCJU-IN 
o 



2 3 4 5 
Population Migratory Score 



Figure ^.ADCYAP1 allele length and migratory status. 

correlation, p > 0.05). 



No correlation was found between allele length and migratory score (Spearman's 




SD(1) LM(2) BC(3) Ml (4) IN (5) MS (6) 



Oregon 



Slate-colored 



Figure 2. ADCYAP1 allele length within sub-species groups. ADCYAP1 allele length did not differ between populations within the 
sub-species for which we have multiple sampling locations. Shown are means and standard errors for Oregon juncos listed by breeding 
grounds (panel A) and slate-colored juncos listed by wintering grounds (panel B). Population migration scores are shown in parentheses with 
higher numbers indicating populations that migrate longer distances (see Table 1 for more details). 
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from Guadalupe island) contained between two and four alleles 
(Table 4). A single allele (276 bp) was the most common in all 
populations except /. p. alticola from Guatemala where the 279 bp 
allele was most common. All but one of the additional alleles (pre- 
sent in a single individual) differed from the most common allele 
by multiples of three base pairs. All populations, both separately 
and combined, were in Hardy- Weinberg equilibrium (all p > 0.05). 



Population mean allele length of CLOCK did not correlate with 
migratory status (Spearman's correlation, rho = -0.38, p = 0.159; 
Figure 3) when all populations were compared. Among the sampled 
populations of Oregon juncos (/. h. oreganus and/, h. thurberi) and 
slate-colored juncos (/. h. hyemalis) there was a significant effect 
of population within subspecies on CLOCK allele length, but not 
a difference in CLOCK allele length between populations (nested 



Table 4. Distribution of CLOCK allele length within each population. Table lists 
the number of individuals, from each population, possessing each of the identified 
CLOCK alleles. Number of individuals genotyped is indicated as "n" for each 
population. Het = calculated heterozygosity for the population. Other abbreviations 
follow from Table 1 . 



Population 


267 


269 


270 


273 


276 


279 


285 


n 


Het 


SCJU-MS 


0 


0 


3 


2 


33 


2 


0 


20 


0.316 


SCJU-IN 


0 


0 


3 


0 


36 


1 


0 


20 


0.188 


SCJU-MI 


1 


0 


8 


2 


49 


0 


0 


30 


0.319 


CRJU 


0 


0 


9 


0 


36 


0 


1 


23 


0.357 


WWJU 


0 


1 


20 


0 


35 


0 


0 


28 


0.490 


ORJU-BC 


0 


0 


2 


0 


27 


1 


0 


15 


0.191 


ORJU-LM 


1 


0 


5 


4 


56 


0 


0 


33 


0.257 


ORJU-SD 


0 


0 


17 


0 


57 


0 


0 


37 


0.366 


GHJU 


0 


0 


3 


1 


52 


0 


0 


28 
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Figure 3. CLOCK allele length and migratory status. No correlation was found between population mean CLOCK allele length and migratory 
status (Spearman's correlation, p > 0.05). 
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ANOVA, subspecies difference: F^^^^ = 1.235, p = 0.267; popula- 
tion within subspecies: F^^^^ = 2.468, p = 0.045; Figure 4). In both 
groups, populations that migrated longer distances possessed long- 
er alleles at the CLOCK locus. 

Individual variation in migratory restlessness 

We tested for an association between mean ADCYAPl allele length 
and level of migratory restlessness within the two captive popula- 
tions of Oregon juncos (/. h. thurberi, Laguna Mountain and UC- 
San Diego). The most informative linear model described indi- 
vidual migratory restlessness in relation to population, mean AD- 
CYAPl allele length, and the population by allele length interaction 
(with Laguna Mountain as reference: effect of allele length = 712.9, 
p = 0.029; effect of population = 117,459.8, p = 0.056; 
population^allele length effect = -732.3, p = 0.053; model = 0.225, 
p = 0.036). This interaction model performed better than all other 
tested models (Akaike's Information Criterion = 669.3; all other 
models AIC > 670.8; see Methods for list of models). Because the 
interaction term indicates a different effect of allele length on migra- 
tory behavior in each population, we ran separate linear mod- 
els for each populations^' which confirmed that longerindividual 
mean ADCYAPl allele length marginally predicted greater migrato- 
ry restlessness in the migratory Laguna Mountain population 
(effect of allele length = 712.9, = 0.169, p = 0.090; Figure 5, 
Supplementary Table 1), but not in the sedentary UC-San Diego popu- 
lation (effect of allele length = -19.36, R^ = 0.001, p = 0.844; Figure 5). 

There was no relationship between mean CLOCK allele length and 
individual migratory restlessness behavior (Figure 6, Supplementary 
Table 1). The most informative linear model described only a 



population difference: the Laguna Mountain population displayed 
more restlessness than the population from UC-San Diego (with 
Laguna Mountain as reference: effect of population = -1299.4, p = 
0.054; R^ = 0. 102). For a more thorough treatment of this population 
difference, see Atwell et al. (in review). This population-only model 
performed better than all other tested models (Akaike's Information 
Criterion = 670.8; AIC for all other models > 671.32). 

Discussion 

Previous studies have related migratory behavior in birds and fish 
to allelic variation in CLOCK and ADCZ4P7. This study examined 
whether variation in these genes also relates to migratory behavior 
in the avian genus Junco. At the level of populations, we found 
that longer CLOCK alleles were associated with longer migra- 
tory distance within two sub-specific groups, /. h. oreganus and 
/. h. hyemalis. CLOCK was not, however, significantly associated 
with migratory distance across the genus as a whole or with indi- 
vidual intensity of migratory restlessness. In the case of ADCYAPl, 
longer alleles were not associated with migratory distance among 
populations within subspecies groups or across the genus. How- 
ever, longer alleles were associated with higher levels of individual 
migratory restlessness in one of two populations studied. Together 
these findings suggest that allelic variation in both CLOCK and 
ADCYAPl may contribute to differences in migratory behavior at 
some levels of analysis (within- subspecies or within-populations) 
but not reliably across the genus. 

The lack of consistent relationships between molecular genetic 
variation and migratory variation suggested by our data are overall 
consistent with previous research on both CLOCK and ADCYAPl, 
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Figure 4. CLOC/C allele length within sub-species groups. Within each subspecies group for which we have multiple sampling locations, 
there is a trend for populations that migrate a longer distance to have longer CLOCK alleles. Figure shows mean and standard error for 
Oregon juncos listed by breeding grounds (panel A; ANOVA, p = 0.056) and slate-colored juncos listed by wintering grounds (panel B; 
ANOVA, p = 0. 1 8; see text for details). Population migration scores are shown in parentheses with higher numbers indicating populations that 
migrate longer distances (see Table 1 for more details). 
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Mean ADCYAP1 Length (bp) 



Figure 5. ADCYAP1 and individual variation. Plot shows mean allele length of ADCYAP1 against individual migratory restlessness score 
for Oregon juncos from the migratory population at Laguna Mountain, CA, USA (solid circles, solid line) and the sedentary population in San 
Diego, CA, USA (open circles, dashed line). Mean allele length marginally predicts migratory behavior in the Laguna Mountain, but not San 
Diego, populations. 



a: 
o 



o 
o 
o 

00 



o 
o 
o 

CD 




272 273 274 275 

Mean length of CLOCK Alleles 



276 



Figure 6. CLOCK and individual variation. Plot shows individual migratory restlessness score for Oregon juncos against the mean length 
of the CLOCK alleles from the migratory population at Laguna Mountain, CA, USA (solid circles) and the sedentary population in San Diego, 
CA, USA (open circles). Points are slightly jittered for each population to improve visualization. There is no significant relationship between 
CLOCK length and migratory restlessness. 



as these genes are related to migratory and breeding behavior in 
only some of the species investigated. Within birds, CLOCK is 
associated with migratory and breeding behavior in blue tits^'^, but 
not black caps, bluethroats, great tits, or swallows Similarly, 
a relationship between ADCYAPl and migratory behavior has been 
found only in blackcaps-^, but has yet to be reported in other species. 
Despite these conflicting results, allelic variation in each of these 
genes is potentially consistent with a direct relationship with migra- 



tory behavior (see below), suggesting that these genes may indeed 
be important to the evolution of migratory and breeding behavior. 

There are at least three reasons why the associations reported here 
might differ according to which populations or species were investi- 
gated: 1) variation in other, related genes (i.e., genetic background); 
2) variation in the degree of linkage between the allelic variation we 
measured and functional genetic differences; or 3) variation in the 
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environment under which the association was sought. These expla- 
nations have been suggested for other genes with variable patterns 
of association with phenotype (e.g., Korsten et al?-). We consider 
how each of these three explanations for variation in relationship 
applies to migratory behavior and CLOCK or ADCYAPl, and we 
suggest that population differences in genetic background are the 
most likely explanation. If this is indeed the case, then future stud- 
ies of candidate genes will benefit from expanding the number of 
target genes examined to include related genes that may contribute 
to differences in genetic background. 

Mechanisms relating allelic variation to migratory phenotype 

CLOCK alleles differ by multiples of three base pairs (Table 4), con- 
sistent with previous results in other avian species (e.g., Johnsen 
et al. Steinmeyer et al^^) and with the location of the repeat with- 
in the coding region of CLOCK^^. It has previously been suggested 
that presence in the coding region might make these polymor- 
phisms functionary particularly because the repeat is responsible 
for a polyglutamine tract expansion in a region known to control 
the rate at which the CLOCK protein activates the downstream cir- 
cadian pathway^'". Polyglutamine expansions control the activity of 
many genes (e.g.. Chamberlain et aV'^). Thus, variation in allele 
length could play a direct role in modifying migratory behavior 
by altering the threshold day length signal to activate or suspend 
migratory behavior, or by modifying its downstream effects. 

The ADCYAPl alleles varied primarily by 2 base pairs (Table 3), 
consistent with previous results and the fact that this microsatellite 
polymorphism falls within the 3' untranslated region (UTR) of the 
ADCYAPl gene^^'^l The mechanism by which this variation may 
relate to behavioral differences remains unknown. Di-nucleotide 
repeats in the regulatory regions of genes, such as 3' UTR, can mod- 
ify the expression of a gene^^. This polymorphism may simply be 
linked with another causative mutation in the coding region-^. The 
amino acid code of ADCYAPl, however, is 97% conserved from 
humans to chickens^^ making 3' UTR based variation in expression 
far more likely than a polymorphic variation in protein sequence. 
Studies of in vivo expression of various alleles in multiple tissues 
would be necessary to confirm a role for this mutation in modifying 
the expression of ADCYAPl in living animals. 

Role of variation in genetic background 

Migratory behavior is a complex phenotype consisting of changes to 
movement, orientation, feeding, fattening, and sleep patterns''\ and 
therefore may be modulated by many genes. The fact that CLOCK 
and ADCYAPl are related to migration in several systems suggests 
that they may be some of the key regulators. However, the large com- 
ponent of migratory behavior that is not explained by these genes 
may be related to other genes in related pathways. Changes in these 
other, related genes may thus reduce or modify the effect of CLOCK 
and ADCYAPl on migratory behavior, and explain the pattern of 
relationships identified in both previous studies ^^'^^'^^'^'^ (described 
above) and this work. Similar patterns have been found for other 
complex traits such as response to stress, including in experimental 
evolution of yeast^'\ Despite similar repeated phenotypes, indepen- 
dently evolved lineages did not duplicate the same transcriptional 
(and therefore genetic) responses to a novel stressor^'l This suggests 
that divergent Junco populations may have utilized unique genetic 



mechanisms to achieve the same phenotypic end, thus masking the 
role of any single candidate gene across the genus. 

Differences in genetic background may explain the lack of genus- 
wide patterns of relationship between migratory phenotype and 
CLOCK or ADCYAPl and the pattern of relationship found in pre- 
vious studies. If other, as yet unidentified, genes also contribute to 
divergence in migratory and breeding behavior, then the relation- 
ships identified between phenotype and genotype would be limited 
to the groups or species in which CLOCK or ADCYAPl are respon- 
sible for the divergence. In contrast, those species and populations 
that do not demonstrate a relationship between migratory pheno- 
type and CLOCK or ADCYAPl may have diverged due to changes 
in other genes. These genes might be involved in the circadian 
pathway, such as BMALl or Period 2^^; the sensing of photoperi- 
od changes, such as thyroid hormone deiodinases^'^ or vertebrate 
ancient opsin^'^'; or transcriptional regulation of these systems^'^. 

Differences in genetic background may also explain the details of 
the relationships found between migratory phenotype and CLOCK 
ox ADCYAPL The Oregon and slate-colored subspecies groups both 
demonstrate a positive relationship between CLOCK allele length 
and migratory behavior. Importantly, the population mean CLOCK 
allele lengths are not different between the subspecies groups, 
and all of the slate-colored juncos migrate longer distances than 
the Oregon juncos. Similarly, allele length of ADCYAPl predicts 
migratory restlessness in the altitudinal migrant Laguna Mountain 
population, but not in the sedentary San Diego population despite 
the fact that both populations exhibit variation at the locus and in 
migratory restlessness. While sedentary populations are known to 
exhibit seasonal migratory restlessness^^ it is possible that changes 
in other genes inhibit migration in this sedentary population suf- 
ficiently to overwhelm the variation in restlessness contributed by 
ADCYAPL Together, these results suggest that variation in other 
genes across the genus account for some of the shifts in migratory 
phenotype, but that CLOCK and may still play important 

roles in further modification of migratory behavior. 

The role of linkage 

An alternative explanation for the finding that allele length of both 
CLOCK and ADCYAPl were related to migratory behavior at only 
some levels (within subspecies or population) of analysis is that 
the alleles we assessed are not functional, but are instead geneti- 
cally linked to functional differences. If so, this linkage may be 
distorted in some populations due to recombination or mutation. 
For those populations, the alleles we measured would no longer be 
informative as indicators of migratory behavior. The presence of 
the relationship in multiple independently derived lineages could 
be explained by independent events that link the allele at the causal 
locus to the same microsatellite alleles. For example, in great tits 
(Parus major), only one of four populations studied, demonstrated 
a relationship between a single nucleotide polymorphism (SNP) in 
dopamine receptor 4 (DRD4) and exploratory behavior^-, perhaps 
due to the breaking of a genetic linkage^''^. 

In the case of CLOCK and ADCYAPl, however, several pieces of 
evidence point against disrupted linkage as the source of variation in 
the genotype-phenotype relationship. Similarly to previous studies 



Page 10 of 16 



FIOOOResearch 2013, 2:1 15 Last updated: 29 JAN 2014 



that found a significant relationship between CLOCK or ADCYAPl 
and migratory behavior (described above), the alleles associated 
with increased migratory behavior in this study were always longer. 
If linkage were responsible for the association, there is no reason to 
predict that longer alleles would repeatedly be linked to increased 
migratory propensity. This does not, by itself, rule out linkage-it 
is possible that this directional pattern is coincidental and will be 
reversed in another species yet to be studied. In addition, unlike 
the length polymorphisms studied in CLOCK, the SNP in DRD4 
related to exploratory behavior is a synonymous mutation^Mhat is, 
it causes no change in the protein that is encoded by the gene. Simi- 
larly, the variation in ADCYAPl is predicted to affect expressions^ 
while the SNP in DRD4 is unlikely to do so (but see Duan et aV^^). 
Given the existence of a known mechanism that could relate allelic 
variation in CLOCK and ADCYAPl to functional consequences 
(detailed above), disrupted linkage seems less likely to account for 
variable strength of association than other possible causes. 

The role of environmental factors 

A final potential explanation for the lack of a consistent relationship 
between allelic variation and migratory phenotype is that the roles 
of CLOCK and ADCYAPl are modulated by the different environ- 
ments that populations and subspecies experience. In three- spine 
sticklebacks {Gasterosteus aculeatus), for example, environmen- 
tal variability in the presence of predators masks genetic variation 
related to several behavioral traits^ ^ Similarly, many environmen- 
tal cues, such as temperature and precipitation, modulate migratory 
behavior'^; such environmental modulators may account for some 
of the differences between junco populations and may mask the 
effects of variation in CLOCK and ADCYAPL 

Environmental differences do not, however, appear to be consist- 
ent with the pattern of relationships found between migratory 
phenotype and CLOCK or ADCYAPl. First, the positive relation- 
ship between individual variation in migratory restlessness and 
length of ADCYAPl was found in only one of two tested popu- 
lations (Laguna Mountain), despite the fact that they were tested 
in an identical common garden experiments^ and winter in the 
same environment^^ Second, the within subspecies pattern seen 
between CLOCK allele length and migratory distance is seen in both 
Oregon and slate-colored juncos, and these two subspecies differ 
in the amount of environmental divergence between populations. 
The sampled Oregon populations breed in divergent climates and 
winter in similar environments, while the slate-colored popula- 
tions all breed in a similar environments, but migrate different dis- 
tances to divergent wintering grounds^^^^i Together, the findings 
in CLOCK and ADCYAPl suggest that environmental differences 
between populations, while potentially important, are not sufficient to 
explain the pattern of associations that we found. 

Conclusions 

We found no evidence for a predictable relationship between 
migratory behavior and the lengths of CLOCK or ADCYAPl 
alleles across the genus Junco. However, we found some support 
for a role of CLOCK allele length within subspecies groups, and 
ADCYAPl allele length within a population, in modifying migratory 



behavior. The lack of consistent detectable associations between 
allelic variation and behavioral variation among populations and 
individuals as reported in our study, adds an important layer in fur- 
ther understanding the potential roles of candidate genes in com- 
plex ecological and evolutionary processes. To explore this layer 
further, future sequencing projects may wish to similarly consider 
both between-group and within-group variability in genotype- 
phenotype relationships when conducting analyses of other candidate 
genes. 

Focusing on variation in a small number of genes with predicted 
functional significance holds the potential to identify major func- 
tional hubs in complex phenotypes across the animal kingdom. The 
redundancy of the genome and complexity of functional genetic, 
regulatory, and signaling networks, however, means that over evo- 
lutionary time scales, candidate loci may gain or lose functional 
significance in the presence of an alternative environment, linkage, 
or genetic background. Our results highlight the fact that differ- 
ent mechanisms may contribute to variation in complex phenotypes 
when examining populations or species. This does not mean that 
candidate gene approaches should be abandoned altogether. Rather, 
with the increasing availability of next-generation sequencing, it 
may become more efficient to expand the scope of candidate gene 
sequencing projects to include other genes in the same pathway, 
as this breadth may provide insight into other, as yet unidentified, 
genetic mechanisms. 
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Supplementary material 



Supplementary Table 1. Individual variation in migratory restlessness and allele length. Allele calls for both CLOCK and 
ADCYAP1 are provided for each individual. Migratory restlessness score provides a cumulative score of migratory behavior 
over the sampling period (March 3 to August 1). See Methods for details. 
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Peterson and co-authors have investigated the association between genotypes of two candidate genes ( 
Clock and ADCYAP1) and migratory behaviour in the avian genus Junco. They examined this association 
among subspecies, populations and individuals. Previous studies have found contradicting results 
regarding the association between phenotypic variation and genotype of these genes; therefore, more 
similar studies are required in order to determine the generality of this relationship and the role of these 
genes. 

The authors have sampled and analysed a large number of individuals and populations and analysed the 
data adequately. Unfortunately, the results are mostly negative and prevent from reaching clear-cut 
conclusions regarding Clock or ADCYAP1 and migratory propensity. However, studying the relationship 
between genotypes and phenotypes on several levels (subspecies, population, and individuals) enabled 
performing interesting comparisons and detecting some positive patterns. 

The authors have properly discussed the results and potential explanations and mechanisms for the 
differences found among the levels examined. 

I agree with the authors that correlating genomes with phenotypic variation may be the next step to 
identify potential candidate genes that may have a role in complicated behaviours such as migration. 

Despite the mostly negative results, I believe that this nicely-written paper provides a fine contribution to 
the literature. One technical comment - it would be interesting to learn the nature of the ADCYAP1 alleles 
that are just one base-pair apart. 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard. 
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Peterson and coworkers investigate phenotypic correlates of genetic variation in two species of the bird 
genus Junco at the level of individuals, populations and species across taxa, which is a valuable 
contribution to our understanding of the genetic basis of the complex behavioural phenotype of migration. 
The methodological approach the authors have chosen is a candidate gene approach at two loci (Clock 
and ADCYAP1), which have previously been analysed in the context of their putative involvement in 
shaping migratory phenotypes in different taxa. The study is a direct follow up from the study by Muller et 
al. 201 0, and a very welcomed attempt to repeat previous findings in another system to obtain a better 
understanding of associations between genetic variation and behaviour, as well as to assess their 
generality. This is much needed since associations between polymorphisms in candidate genes and 
behavioural traits have often proven to be difficult to repeat across populations or species, and replication 
is indeed a valuable tool to infer the generality of association patterns. 

The difficulty and my primary concern with the study as well as the precursor study is that candidate 
genes (defined as genes that have been studied in one (model) organisms and there have been found to 
influence the expression of one focal phenotype) for migratory traits are probably not know to date, as 
none of the model organisms show a migratory phenotype. The genes in focus here are candidate genes 
for candidate traits, probably involved in shaping the migratory phenotype (i.e. the circadian phenotype 
involved in timing of migration). Thus, conclusions drawn upon these data based on correlative evidence 
should be evaluated with great caution. 

This study is useful in the way that it adds data to the pool of studies that have explored candidate genes 
of candidate traits for the migratory phenotype, which are difficult to interpret and don't allow drawing any 
clear conclusions upon. The study is well presented and the additional data confirm the apparent 
inconsistency and the difficulty to interpret these data, probably highlighting the weakness of the 
candidate gene approach in this field of migratory research, and I support their publication. With 
candidate gene studies per se we do not expect to get huge effects, and in the area of migration genetics, 
I suggest that future research should avoid using a candidate gene approach, as this is unlikely to yield 
clear results - at least unless clear candidate genes for the migratory phenotype have been identified de 
novo. 

I have read this submission. I believe that I have an appropriate level of expertise to confirm that 
it is of an acceptable scientific standard. 
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